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1 Introduction 

Quantum field theories admit a variety of operators defined not by insertions of the fundamental fields, 
but by constraints which change the domain of the path integral in field space. An operator defined 
by such a prescription is called a defect operator. A famous example is the twist operators of 2d 
conformal field theory. More generally, one can define a defect operator inserted along a submanifold 
L by deleting L and requiring the fields to have prescribed singularities as one approaches L. The 
effect of the insertion can be "measured" by evaluating the expectation values of ordinary operators. 
When one can detect the presence of the defect from afar (for example, because some field strengths 
are now required to belong to a nontrivial cohomology class), the insertion is said to have created 
topological disorder. 

The first example of a defect operator in gauge theory is probably the 't Hooft loop operator in 
4d gauge theories [1] which can be used as an order parameter for Higgs phases. It also plays an 
important role in the context of electric-magnetic duality of J\f = 4 gauge theories in 3 + 1 dimensions 
[2-4] . The duality exchanges states with electric and magnetic charge and therefore exchanges Wilson 
loop operators and 't Hooft loop operators. 

There is a somewhat similar story for 3d gauge theories. In such theories there is a often a duality 
which exchanges elementary excitations in the Coulomb phase with Abrikosov-Nielsen-Olesen vortices 
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in the Higgs phase. A 3d operator creating a very heavy vortex with world line L is a defect loop 
operator and may be regarded as analogous to the 't Hooft loop operator in 4d creating a monopole. 
Such a defect operator is defined by the fact that the gauge field has a fixed holonomy around any 
small loop linking L. In 4d an analogous construction gives a surface defect [5]. 

It should be noted that the definition of the vortex loop operator is independent of the existence of 
the Higgs phase, or vortex solutions, or even of a dynamical gauge field. The definition makes perfect 
sense in the topological pure Chern-Simons theory as a defect in the gauge connection, although the 
defect, in that case, can be identified with a Wilson loop [6, 7]. In a theory possessing an abelian 
global symmetry, a vortex loop can be defined by gauging a global symmetry using a non-dynamical 
flat connection satisfying the holonomy condition. It is also not necessary that the defect be defined 
on a closed loop, however, to preserve gauge invariance, an open contour must extent to the boundary 
of spacetime. 

Much more can be said about defect operators when the theory is supersymmetric. All of the 
operators mentioned above have BPS analogues in supersymmetric theories in 2 + 1 and in 3 + 1 
dimensions. Extending the definition of a defect so as to preserve a fraction of the supersymmetry 
can require imposing conditions on additional fields. This is analogous to the inclusion of fields other 
than the connection in the definition of a supersymmetric Wilson loop. In this work, we will define 
the supersymmetric analogue of the vortex loop. Exact computation of the expectation value for a 
supersymmetric defect may be feasible by employing localization techniques ([8-10]). This was carried 
out for the supersymmetric version of the 't Hooft loop on S 4 in [11]. Here, we extend previous results 
for localization of supersymmetric gauge theories in 2+1 dimensions (see [12] for the original derivation 
and [13] for a review) to include the supersymmetric vortex loop. We will also discuss the role played 
by supersymmetric vortex loops in the context of mirror symmetry. 

In section 2, we define several versions of the abelian vortex loop. We discuss the operator's trans- 
formation under Witten's SL(2, Z) action on conformal field theories. We then extend the definition to 
accommodate supersymmetry. In section 3, we employ localization to evaluate the expectation value 
of the supersymmetric vortex loop for a generic superconformal Af = 2 gauge theory on S 3 . The result 
can be inferred from the SX(2, Z) action. We provide an independent derivation using the original 
definition. In section 4, we demonstrate, with a few examples, the role of the vortex loop in IR duality 
of Af = 2 gauge theories. We conclude with a discussion of possible extensions. 

2 Definition of the abelian vortex loop 
2.1 Vortex loop in a gauge theory 

Defect loop operators in 3d gauge theories have been previously introduced in the context of Chern- 
Simons theories [6, 7[. Such a loop operator is specified by giving a loop 7 in the 3-manifold M and 
an element /3 of the Lie algebra g of the gauge group G. The holonomy around any small loop linking 
7 of the gauge connection A is required to approach /? as we shrink the loop size to zero. With this 
condition, A is singular on 7. Separating A into a smooth part A' and a singular part A", we can 
write 

F A „=j3*fr] (2.1) 

where -k[y] is a 2-form current supported on 7 whose cohomology class is the Poincar'i^ce dual of 
[7] G H\(M,Z) (we assume that M is orientable). In [5], the authors defined surface operators in 
Af = 4 SYM theory in 3 + 1 dimensions using a similar prescription. The prescription for a subset of 
these codimension 2 operators, the one in which only the connection is singular, coincides with the 
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definition above (substitute a for /3). As noted there, the data specifying the singularity is actually 
only e I/3 , and equation 2.1 should be handled with care. If the gauge group G is U(l), ft is simply a 
real number. To avoid ambiguity it is sufficient to restrict the range of f3 to the interval (— 7r,7r). In 
what follows, we will sometimes assume this restriction. We will also set 

? = 2ttA q€ (-l,Vj (2.2) 

From now on we will assume that G is abelian and will call such a loop operator a gauge vortex 
loop. 

2.2 Global vortex loops and the SL(2, 1) action 

We can define a similar loop operator if G is a global symmetry group rather than a gauge symmetry 
group. To this end we merely set the smooth part A' of the gauge field to zero and set the singular 
part (A" above) to be a fixed flat connection whose holonomy along a loop linking 7 is e 1 ^ . We will 
call such a loop operator a global vortex loop. It appears naturally when we consider the action of 
Witten's SL(2, Z) on Wilson loops. 

Following [14], we consider a conformal field theory in 2 + 1 dimension with a choice of an abelian 
global symmetry current J. We couple J to a background gauge field A, and consider the partition 
function as a functional of A 

Z[A}= j 2 3 $ e ^[*]+ l / v ^<i 3 x,/^ tl +... ( 23 ) 

where refers to seagull terms necessary to ensure invariance under gauge transformations of A. 
Here <E> and S[&] are short hand for the fields and action of the theory. In addition to the seagull 
terms, one can add extra terms which are gauge-invariant functionals of A alone. If we wish to preserve 
conformal symmetry, a natural choice is the abelian Chern-Simons term 
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AAdA (2.4) 



If we want to associate A to a non-trivial principal bundle, then this is only defined for integer a. 
However, for now, we will assume that A is a connection on a trivial bundle, and allow arbitrary real 
values of a. To the triplet (S[$], J, a) we associate the following functional of A: 

Zj a [A] = J J)^^^^! Va^xJ^A^ + .^ + i^ J AAdA 

Witten defined an SL(2,Z) action on such triplets [14]. The action of the T generator is merely 
a shift of a by 1. The action of S is defined by first promoting A to a dynamical gauge field and then 
replacing J with the topological U(l) current for this new gauge field 

J top =—-kdA (2.6) 

Together, these operations generate an action of SL(2, Z) on the set of triples. At the level of partition 
functions, these operations act as follows: 

(T-Zj, a )[A] = Zj, a+1 [A] (2.7) 
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(S ■ Z Ja )[A] = I T>A'T>$e iS ^ ]+i J ' Vad 3 xj^A' )l +...+^ j ■ A'/\dA'+^ J • aa&a' ( 2 .8) 

Our first goal is to extend these operations to observables more general than the partition function 
Zj.a- Specifically, we would like to consider insertions of Wilson loops for the background gauge field 
A. First, it is useful to define a slight generalization of an abelian Wilson loop. Recall that the abelian 
Wilson loop operator is specified by a loop 7 : S 1 — > M and a charge q € R and is defined as an 
insertion of 



e ig ^ A (2.9) 
into the path integral. This can be rewritten as 

e^^ uAA (2.10) 

where a; is a closed 2- form current with support along 7, defined so that the integrals in the previous 
two expressions agree for all 1-forms A. In local coordinates (r, 8, z) where the loop lies along the z 
axis, we can write: 



u = 2Tiq8 1 = qS(r)dr A dO (2.11) 

Such a term can simply be added to the action, since in the abelian case we do not have to worry 
about path ordering the exponential. 

We can consider such an insertion for more general 2-form currents lj. Invariance with respect 
to infinitesimal gauge transformations forces ui to be closed. Invariance with respect to "large" gauge 
transformations (i.e. gauge transformations which are topologically nontrivial maps from the abelian 
Lie group G to M) requires the de Rham cohomology class of uj/2tt to be integral. This means that 
to arises as the field strength of some connection A u on a U(l) bundle, and we can write the Wilson 
loop as a BF coupling to this new background gauge field: 

In the case 

w = 2%q5 1 (2.13) 

where 8 1 is a 2-form current supported on 7 which is Poincar'i^ce dual to 7, the integrality condition 
on oj reduces to the requirement that the class ^[7] £ Hi(M,M) is integral. In particular, if 7 is 
homologically trivial, there are no integrality constraints on q. Near 7 the connection A u in suitable 
coordinates looks as follows: 

A, = qdO (2.14) 

More generally, if we take ui to be supported in a small tubular neighborhood of 7, we get a regular- 
ization of the Wilson loop along 7. 

With this in mind, we define a insertion of W w by 

(W u ■ Z. la )[A] := J V ^Sm^.[ V9d 3 ^A, + ... + ^ f AAdA+^ f A^AdA 

Let us see how the SX(2,Z) generators act on it. It is clear that T commutes with W u , so we only 
need to consider 
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({Sr 1 W u S)-Z Jta )[A] = J T>AiVA 2 V$exp(iS[$]+i J ^d 3 xJf"A lt , + ... + ^-J AiAdAi+ (2. 



16) 



+ h I A2 A dAi + i J Au A dA2 _ h ] AA dA y (2 ' 17) 

We can see A2 enters only via a term J (A\ — A u — A) A dAi , and so the integral over A2 produces 
exactly a delta function setting A\ = A u + A [14]. This leaves 

j V<$> exp [iSW\ + * f y/gdPxJUAup + A„) + ...+ '£ J (A + A u ) A d(A + A^ (2.18) 

= Z Jta [A + Au] (2.19) 
Note that, even at A = 0, this gives an insertion of an operator: 

cxp (i J s /gd 3 xJ"A Ufi + ... + H y A„AwJ (2.20) 

In the case where u = 27r<j(5 7 , A m is a flat gauge field with a holonomy e 27 ™ 9 , and this operation has 
the same effect on the path integral as prescribing that all fields charged under the current J pick up 
a fixed monodromy around the loop 7. There is also an additional a-dependent phase factor related 
to the self- linking number of the loop (which must be regularized by specifying a framing). We see 
that in this special case the operation S maps the charge-g Wilson loop for the background gauge field 
A to the global vortex loop with holonomy e 27 ™ 9 . More generally, for arbitrary w satisfying the above 
integrality conditions we can define an operation 

(D u -Zj ia [A]) = Z Jta [A + A u) ] (2.21) 

What we have demonstrated is that the operation S maps W u to D w : 

S^WuS = D u (2.22) 

More precisely, the equality holds up to a phase factor which depends not only on u>, but also on the 
Chern-Simons coupling of the theory on which these operations act. 

One can rephrase this result in terms of Wilson loops for dynamical gauge fields. Gauging a 
symmetry (without adding a Chern-Simons term) is the same as applying the operation S. The 
resulting theory has a new global symmetry U{l)j whose current is -kdB/2ir, where we denoted by B 
the dynamical gauge field, to distinguish it from the background gauge field A which couples to the 
U{l)j current. Noting that S 2 = C (the charge conjugation), we get 

SW^ = D- U S. (2.23) 

Applying this to the partition function z U j"^ a9ed [B\, we learn that 

< W w [B] >= • Z gauge [A] = Z gauge [A-AJ\. (2.24) 
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In particular, setting ui /2?r to be the delta- function supported on a loop 7, we see that a global vortex 
loop for the U(l)j symmetry is nothing more than an ordinary Wilson loop in the underlying gauge 
field. Similarly, we find 

< D U [B] >= W u ■ Z„[A] (2.25) 

which shows that, in the absence of a Chern-Simons term, the gauge vortex loop by itself is somewhat 
trivial: it merely modifies the functional dependence of Z gauge [A] on the background gauge field A 
which couples to the U(l)j current. 

2.3 Pure Chern-Simons theory 

Before moving on to the supersymmetric version of the vortex loop, let us briefly comment on how 
the gauge vortex loop behaves in pure bosonic Chern-Simons theory. It was argued in [6] that such 
a defect operator should be equivalent to a Wilson loop. However, we have seen above that for an 
abelian gauge group the gauge vortex loop is somewhat trivial, its only effect being a modification 
of the U(l)j current by a c-number term supported on the loop. To see that this agrees with the 
behavior of the Wilson loops, we recall the formula for the expectation value of a product of Wilson 
loops in U(l) Chern-Simons theory at level k [7] 

(Y[exp(iq a [ A)) =exp(^Y,1^(7a,7b)) (2.26) 
•'To V K a,b 7 

where 3>(7 a ,7&) is the linking number of the loops a and b. The latter can be written in terms of the 
corresponding gauge fields A a , with dA a — uj a the 2-form delta function supported on j a , as 

$(7a,7&) = J A a = ^ J A a Auj b (2.27) 



Letting 



a 

we can rewrite the expectation value as follows: 



(Y[exp(iq a f A)) = cxp (J^ f A * A dA u J . (2.29) 

Thus the insertion of a collection of Wilson loops is equivalent to a phase factor which depends on k 
as well as on a flat connection A u . This is compatible with the claim that a collection of Wilson loops 
in pure abelian Chern-Simons theory is equivalent to a gauge vortex loop D^, for some 2-form current 
u>, which in turn is trivial up to a phase. 

2.4 Supersymmetric vortex loops 

We would like to extend the considerations above to the supersymmetric case. Specifically, we will 
work with theories with Af = 2 supersymmetry (4 real supercharges). It is convenient to work in 
M = 2 superspace, with fermionic coordinated 9 a and a superspace derivative D a . For the theories of 
interest, the dynamical fields can be organized into chiral and vector superfields. The gauge field is 
part of a vector multiplet, and all fields in this multiplet take values in the adjoint representation of 
the gauge group. A vector superfield V satisfies V = and contains a vector field A„, a real scalar 
<7, a complex spinor A Q , and a real auxiliary scalar D. Matter fields live in chiral multiplets and take 
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values in some representation of the gauge and flavor groups. A chiral superfield $ satisfies D a Q = 0, 
and contains a complex scalar 0, a complex two component spinor ip a , and an auxiliary complex scalar 
F, 

We are interested mainly in a class of renormalizable gauge theories with abelian global symmetries. 
These are defined by a UV action which includes a kinetic term for the matter fields of the form 

^charged matter kinetic = - / d 3 xd 2 8d 2 0j2 (ff^^i) (2-30) 

i 

= 12 f^ x (( D ^) t ( D ^Y + i^D^ + F l F l - &<rV + faDtf - + i^X^ - i^Xcff) 

i •* 

(2.31) 



and a supersymmetric Yang-Mills action for the fields in the gauge multiplet 

1 / JzfOfSTr, Q 



S Y a„g Milk = 4 / d 3 xd 2 6d 2 6Tr f (]z 2 ) (2.32) 



= i j d 3 xTr f (^F^F^ + D^uD^u + D 2 + iXj^D^x) (2.33) 
where £ is a linear multiplet defined by 

£ = D a D a V (2.34) 

£ f = £ (2.35) 

D a D a Y> = D a D a E = (2.36) 

In addition, one can allow supersymmetric completions of Chern-Simons terms. In the abelian case, 
these have a very simple superspace expression 

k f 

S'abelian Cliern Simons = ^~ / d 3 xd 2 0d 2 9Trf (VE) (2.37) 

■ y d 3 xT r/ (^"A^Ap - AA + 2Dct) (2.38) 



A- 
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As a slight generalization of this, we can also consider a "off-diagonal" Chern-Simons term, also known 
as a BF term, coupling two or more different abelian gauge fields. The supersymmetric completion of 
this has the form 

SbF = / d3 ' xd2dd:2 ^ Tr f (E 1 ^) (2.39) 



/ d 6 x \E^ p A\d v A l p - -X ] X l + D l a 3 ) (2.40) 

We will consider the S 3 partition function for these theories with insertions of Wilson loop and defect 
loop operators. The supersymmetry transformations on S 3 are generated by Killing spinors [12]. We 



-7- 



will need the transformations of the vector multiplet fields under the S 3 supersymmetry generated by 
a particular killing spinor e 



M M = -|AV (2.41) 

8a=--\*e (2.42) 

SD = -|iVVe + ^[A f , a]e + \x^e (2.43) 

SX = (- \e^ p F^ lp -D + »7"D m ct - a)e (2.44) 
(5A f = (2.45) 

The fermionic symmetry generated by e will be used in 3 to compute the expectation value of the 
defect operator on S 3 using localization. 

We now attempt to generalize the SL(2, Z) action to the supersymmetric case. To start, it is 
natural to define T by simply adding a supersymmetric Chern-Simons term instead of an ordinary 
one. For S, we should use the supersymmetric version of the BF term and integrate over the entire 
background vector multiplet. It is now a simple exercise to check that the SL(2,Z) relations remain 
satisfied for this generalization. We will not use the (ST) 3 = C property, so we omit a check of 
that relation, but it will be important that S 2 = C, so let us sketch the argument. Consider the 
supersymmetric version of 2.16 defined via the action 2.39 where all vector fields have been extended 
to TV = 2 vector multipelts V, V\, V 2 and V u . When we integrate over the second gauge field, the BF 
coupling gives us a delta function constraint imposing that the first gauge field is the negative of the 
background gauge field, as before. In addition, one can see that the integration over the auxiliary fields 
in the second vector multiplet imposes a similar constraint on the auxiliary fields of the first vector 
multiplet. Thus we see that the net effect is to set the first vector multiplet equal to the negative of 
the background vector multiplet. 

Now consider the supersymmetric generalization of the abelian Wilson loop operator. For a certain 
class of loops 7 preserving supersymmetry, this has the form 



(iq(£{A-iad£)) (2.46) 



exp iq 

For example, on S 3 , the loops must be great circles which are fibers of the Hopf fibration. This 
operator is then invariant under the supersymmetry generated by e. We define an operator W 7 
inserting a supersymmetric Wilson loop in a background vector multiplet as before. Finally, we define 
the supersymmetric vortex loop by the prescription _D 7 = S~ 1 W~ ( S '. 

As before, we can integrate out the auxiliary vector multiplets V± and V 2 to obtain a description 
of the defect terms of the original fields alone. To start, let us write the part of the action containing 
V 2 . In terms of component fields the action (in Euclidean signature) looks as follows: 

S[$,V 1 ,V 2 } = ... + ^ J(-A 2 AdA 1 ) + ^- J d 3 x(-a 2 D 1 -a 1 D2 + ^(X\X2 + XlX 1 )) + tq J (A 2 -ia 2 d£) 

(2.47) 

For simplicity we set to zero the background vector multiplet. We see that the integrals over D2 and 
A2 set <7i and Ai to zero, while the integrals over A2 and 02 impose the constraints 
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dAx = 2irq5~ ( , -kD x = -2mq8 1 A di. (2.48) 

Here di is the volume 1-form on 7 and 8 1 is the 2-form Poincari'^ce dual to [7], as before. Note that 
Di is purely imaginary, which violates the usual reality condition on D. 

As before, we can generalize this operator to account for more general background gauge multiplet 
configurations. We to specialize to S 3 , and pick a supercharge S corresponding to the Killing spinor 
e. Then the BPS condition for an abelian vector multiplet is [12] 

= (iy»(-*F li +d lt tr)-(D + o))e (2.49) 

We would like to find configurations for which only F and D are non-zero. Using u Al 7 A1 e = e, where 
is the Killing vector along the Hopf fibration as in [12], we see we can take: 

F = 2nf(x) *v, D = -2mf(x) (2.50) 

for a function / : S 3 ~ > R. The normalization is for later convenience. Note that the Bianchi identity 
implies 

= ^-dF = d(f * v) = df A -kv = ★Kd,,/) (2.51) 

so that / must be constant along the fibers of the Hopf fibration. Equivalently, we impose that / 
arises from a function / : S 2 — > K by: 

/ = /ott (2.52) 

where tt : S 3 — > S 2 is the projection map of the Hopf fibration. Thus, the operator is really labeled 
by the function / on S 2 . When / approaches a delta- function on S 2 , the corresponding operator 
approaches the supersymmetric vortex loop wrapping a fiber of the Hopf fibration. 

If we apply this construction to the U(l)j symmetry, we get an operator which is a natural 
generalization of the supersymmetric Wilson loop: 

Of = exp(ij y/gd^xfix^A^-ia)^ (2.53) 
That this operator is supersymmetric follows from 

d(v^Af, -ia) = (2.54) 

as can be easily checked. The condition v^d^f = is necessary for gauge-invariance. Taking / to be 
a delta function on S 2 , so that / is a delta function supported along a great circle in 5* 3 , one recovers 
the ordinary supersymmetric Wilson loop. Note that we have picked the normalizations so that when 
/ is a delta function on S 2 integrating to 1, so that / integrates to 2ir on S 3 , we recover the charge 1 
Wilson loop. 



3 Localization in the presence of a vortex loop 

In this section we compute the expectation value of the global vortex loop on S 3 by localization. The 
global symmetry is assumed to be a U(l) subgroup of the flavor symmetry group. We will present 
three approaches to the calculation which yield the same result. Applying the localization procedure 
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in the presence of the defect requires some regularization and the agreement of the approaches pre- 
sented below gives us confidence in the validity of the computation. We begin with a quick review of 
localization for 3d gauge theories. Additional details can be found in [12]. 

3.1 Localization of 3d gauge theories 

The expectation value of BPS operators in 3d M = 2 superconformal gauge theories can be computed 
by localization on S . The relevant result, derived in [12], is that deformation invariance allows us 
to reduce the computation of the infinite dimensional path integral with BPS operator insertions 
to a matrix model with the integration domain given by the Lie algebra of the gauge group. The 
data entering the computation is a UV action with gauge group G, Lie algebra q and Chern-Simons 
levels ki , a set of representations Ri for the chiral matter multiplets and the IR conformal dimensions 
(equivalently R-charges) A,; for each chiral multiplet. The integration measure for the matrix model 
is then 

Vo!(G) da|QeAd(0) (3 ' 1} 
The contribution of a level k Chern-Simons term (for a simple gauge group factor associated to a) is 

e -iirkTr(a 2 ) (g 2 ) 

A Fayet-Iliopoulos term with coefficient r\ contributes 

e 27rir?Tr(a) 

Every dynamical gauge multiplet contributes 



(3.3) 



z gau„ lt iplet (a) = detAd(B) (2sinhM) = JT 2smh(7rp(a)) (3.4) 

pSroots(g) 



and every dynamical chiral multiplet contributes 

^chiral multiplet/ an AetO F 



where p are the weights of R and 

t{z) = -zlog (1 - e 2 -) + 1 (ttz 2 + i L i 2 (e 2 -)) - g (3.6) 

z(p(a),A) = ip(a)-A + l (3.7) 

abelian flavor parameters can be incorporated by shifting p — > p+m. The insertion of a supersymmetric 
Wilson loop in a representation R gives a factor of 

W(a) = ^ f) Tr R (e*™ ) (3.8) 

which for an abelian Wilson loop of charge q reduces to 

e 27Tqa (3.9) 

Integration with the measure 3.1 of the product of all relevant contributions yields the exact expecta- 
tion value. 
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In computing 3.5 we have implicitly assumed a standard 5 exact term, and hence standard kinetic 
operators (Of, Ob), for the fields in the chiral multiplet [12, 15]. In the presence of the flavor vortex 
loop, the kinetic term of a charged chiral multiplet is altered by a background gauge field created by 
the loop. The new term and the revised 1-loop contribution are derived below. The contribution of 
the vector multiplet is unaffected because it is not charged under flavor symmetries. 

3.2 Method 1: using the SL(2,Z) definition of D 

The simplest way to extract the effect of inserting a supersymmetric defect line operator is by using 
the definition of the operation D as 

D q = S^WqS (3.10) 

where D q is the vortex loop with holonomy exp(27rig). There are no integrality constraints on q because 
the large circle on S 3 is homologically trivial. Since we can perform the operations on the RHS at the 
level of the matrix model, it should be possible to compute the LHS indirectly this way. Explicitly, 
suppose we compute the partition function as a holomorphic function of a flavor deformation m (and 
possibly other parameters which we suppress): 

Z(m) (3.11) 

This is the analogue of Z[A] in section 2. The operation S then tells us to treat this flavor parameter 
as a gauge parameter, and integrate over it, with a coupling to an FI parameter rj: 

(S ■ Z)( V ) = f dmZ{m)e 2lTl7lm (3.12) 



The operation Wi tells us to insert a charge-g Wilson loop in the background field corresponding to 
77: 

(W q S ■ Z){rj) = e 27r<? " f dmZ(m)e 2mT,m (3.13) 



Finally, S 1 tells us to integrate over 77 and insert a new FI term, which we will denote m', with the 
opposite sign: 

(S _1 WiS ■ Z){m!) = J e -»*W e 3*OT j dmZ{m)e 2 ^ m (3.14) 

Now to integrate out the variables rj and m, we simply note that the integral over r\ imposes a delta 
function which sets m = m! + iq. Thus we are left with: 

(D q ■ Z)(m) = Z(m + iq) (3.15) 

Indeed, this result can be inferred from the 4d perspective by considering Wilson and 't Hooft loops 
ending on a 3d boundary [16, 17]. This argument was rather indirect; it also raises the question about 
the interpretation of poles in the partition function for special values of m + iq. We now proceed to 
present two more explicit derivations of this result. 
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3.3 Method 2: smearing the defect 

We return to the smeared Wilson loop 



Oj = cxp (i J ^gd 3 xf{x){v ti A ll - ia)j (3.16) 

where / is some real function on S 3 constant along the fibers of the Hopf fibration, as above. It is 
convenient to decompose f(x) as 



f(x) = ±+f (x) (3.17) 

where 

f ^gd 3 xf o = (3.18) 
Js 3 

and q is constant, specifically, q = J s3 y/gd 3 xf. Note that this normalization agrees with the case 
where / is a delta function supported on a great circle, since the integral of / should give 2tt times 
the charge of the Wilson loop. 

For general /, we can decompose Of = O q Of o , so it suffices to study them separately. Actually, 
we will find that it is only O q which contributes to the localized path-integral. Specifically, we claim 
that the operator Of o can be absorbed into a shift of the action by a total S- variation. To see this, let 
us pick a function g : S 3 — > R and consider: 



S( I v /5d 8 xeV(M A ) = / s/9d 3 xJj»(d fl g)(i 7 v (-l;e„ p<T FP'> + d v o) - (D + a))e) (3.19) 



^gd 3 x(d„g) (i(g^ + i^v^—^F^ + d v a) -v„{D + a)j (3.20) 

The term involving D + a is proportional to v^d^g and vanishes if we impose that g, like /, is constant 
along the fibers. The remaining terms can be expanded to give: 

j y/g#x(dpg) ( - l -e^F pa + F^v v + id»a - ie^d v av p ^j (3.21) 
Integrating by parts, the first and last terms can be seen to vanish, and the others give 



^d 3 x^(^(v^g - v v d„g)) - io(V 2 g) 



(3.22) 



The quantity multiplying A 11 can be expanded as: 



V%(9„<?) + v p V 2 g - (V« w )5 M fl - v u V u d p g (3.23) 
Using Vufy = £u,vpv p , this can be simplified to: 

^V 2 S -V^t)„Vs) (3.24) 
The second term vanishes when we impose that g is constant along the fibers, and we are left with: 

i/gd 3 x(V 2 g)(v^A p -ia) (3.25) 
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which agrees with the exponent of the operator Of o above, provided we can find a g such that: 



fo = V 2 .g. (3.26) 

This clearly requires J g3 yfgd 3 xf = 0, so that one cannot use this trick to remove the constant part 
of /. However, if this condition is met, then the equation can be solved, and g will indeed be constant 
along the fibers. 1 This proves that the non-constant part of the operator can be discarded, as <5-exact 
terms do not affect the path integral, i.e. we can replace a Wilson loop localized on a loop by one that 
is smeared uniformly over the entire S 3 ; these differ only by a 5-exact insertion. 

Thus without a loss of generality we can restrict to the case / = ^ (a constant). Then the 
background vector multiplet we must couple the flavor symmetry current to is given by: 

F = 2q*v, D = -2iq (3.27) 
Since dv = 2 * v, the corresponding gauge field can be taken to be 

A = qv (3.28) 
The 5-exact gauged action of a chiral multiplet of conformal dimension 1/2 is given by [12] 

S s = J Jgd^xU^-D^ + a 2 + iD + ^)(f) + ^(i 7 »D^ - ia)ip + F*FJ (3.29) 

Let us couple this to an ordinary BPS background, with a = —D = er , as well as the background 
vector multiplet described above, with A = qv and D = —2iq. We find: 

Ss = J V^ 3 :r(y(-V 2 -2i^d M VW-^o+2 g +^)<^ (3.30) 



The bosonic operator has the form: 



O b = -V 2 + aiv^df, + b (3.31) 



with a = —2q and b = q 2 + a 2 — ia Q + 2q + |, which has determinant [12] 



A/detOs = Y[ ( [£(l + 2)~2am + bj ) (3.32) 

£=0 ^m=-£/2 



i+1 



1=0 v m=-£/2 

Meanwhile, for the fermions, the operator has the form: 



e+i 



I ( [l(i + 2) + 4qm + q 2 + a 2 - io + 2q + - ) ) (3.33) 



O f = i7 /i V M + ic7% + d (3.34) 

with c = — q and d = —icr , which has determinant: 

x This can be seen most easily by working with a mode expansion f = m,n^l,m ni where Yt^ m ^ n 

are spherical harmonics on S li , satisfying V 2 Yf jmiTl = — 1(1 + 2)Y^ jm ,„ and n^c^Y^m^ = tmV^ m „ . Then g = 
m n £/(_i_2) c t,m,nYt,m,n> which is well defined since / has no £ = component by assumption, and, like /, has 
c e,m,n = for all m ^ 0. 
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oo , 1/2-1 , \\ e+1 

I ((£+q+ia +3/2)(£-q+i<7 +3/2) (e(l+2)+4qm+2q-ia +q 2 +a 2 +3/4)) (3.36) 

£=0 ^ m=-t/2 ^ ' ' 

We see that most modes cancel, and we are left with: 

dctOf _ yj ( {£ + q + i<j + 3/2){£-q + t<j + 3/2) 



~ V^tO^ ~ fi V W + 2) - 2g£ + 2g - ia + g 2 + a 2 + 3/4) J ^ ' J 



iiu^+i+8+*o-o + |)(^+i+«-^o-|)y 

- n U-^ +^)+i/2j (3 - 39) 

Note that a Q and a appear in the combination a + iq, so that we can just make this replacement in 
the ordinary one-loop determinant to find: 

Zi- loop = e^+^ (3.40) 

In fact, this computation goes through with minimal changes for chiral multiplets of arbitrary dimen- 
sion (one merely shifts a by an imaginary amount), and we find: 

Z^ioo P = e^- A+ ^ (3.41) 

Note we have obtained the same result as in the indirect argument above. We will now turn to an 
even more explicit argument, where we do not smear out the defect over the sphere but instead work 
directly with a (regularized) delta function background. 

3.4 Method 3: explicit computation in a singular background 

Let us now focus on the specific case where the function / is a delta function, corresponding to the 
dual of an ordinary (unsmeared) Wilson loop. Although we have argued that one can replace the 
delta function by a constant which has the same integral over S 2 , we would like to gain better physical 
insight into the vortex loop by explicitly finding the modes in a singular background. For simplicity 
we will focus on the case where the matter has canonical dimension, although it is straightforward to 
generalize this. 

Recall that the smeared vortex loop on S 3 can be obtained by coupling to a background F = 2irf*v 
and D = —2irif, where / is some function on S 3 constant along the fibers of the Hopf fibration. Here 
we compute the modes and 1-loop determinant explicitly in the case where / is a (infinitessimally 
smeared) delta function supported on a single fiber. 

We will work in toroidal coordinates on S* 3 , with 77 £ [0,7r/2] and 8 and <j> in [0, 2tt). Explicitly, 
we can relate these coordinates to the unit sphere S 3 C M 4 via: 



- 14 - 



:r = cos 77 cos y — cos 77 sin 0, z = sin 77 cos cf>, w — sin 77 sin <fi (3.42) 

The surfaces of constant 77 are torii, which degenerate to great circles at 77 = 0, ^. 
The usual round metric takes the following form in these coordinates: 

ds 2 = drf + sin 2 T]d6 2 + cos 2 r/dcj) 2 (3.43) 

rv2 Id d 1 d 2 1 d 2 , nAA , 

= ~ — smyycosT?— + , 2 + — — ? (3.44) 

sm 77 cos 77 077 or] sin 77 or cos z r\ o<p 

The Killing vector v is given in these coordinates by: 



_ d_ d_ 

or, as a 1-form, by: 



(3.45) 



v = sin 2 rjdO + cos 2 rjd(f> (3.46) 

It satisfies: 

dv = 2 cos 77 sin 77^77 A d6 — 2 cos 77 sin 77^77 Ad<f> = 2-kv (3.47) 

We take the defect to be supported on the great circle at 77 = 0. We define a regularized delta function 
supported on the loop by: 

= g( V /e) 
2ite sin 77 cos 77 

where g(x) has support in < x < 1 and integrates to 1 on [0, j-]> so that / integrates to 2ir. This 
approaches a delta function supported on the great circle at 7/ = for e — >• 0. Then the background 
field-strength we need to consider is 

F = 2irq f * v = -g{i 1 /e)dr t A (d6 - #) (3.49) 

e 

which is solved by a background vector potential 

A = qG(r)/e)d0 - q(G(r]/e) - l)d<f> (3.50) 

where G = g with G(0) = 0, so that G(x) — > 1 for large x, and we pick the constants so that this is 
everywhere well-defined. We also have a background auxiliary scalar 

B = -2^ = - * gg(??/e) (3.51) 
e sin 77 cos 77 

3.4.1 Bosons 

First consider the bosonic operator on S 3 : 

O b = -D„Df + <j 2 + iZ?+ I (3.52) 

We will couple to a defect background, as above, as well as an ordinary a = —D = a Q background. 
Using the form of D above, along with: 
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A»d u 



A 2 



q G(r,/e) d q(G( V /e) -1)5 



sin 77 06 cos 2 77 dip 

q 2 G( V /e) 2 , q 2 (G( V /e)-l) 2 



sin 2 77 



cos 2 77 



(3.53) 



(3.54) 



and looking for an eigenfunction Ob4> = A0 of the form cf> = f^ e tme + m 4> ; we g e ^ the following 
equation: 

1 d / . df\ f(m + qGf (n-q(G-l)) 2 qgMe) 3 , . \ , n 

— : 3- sinTycosTy-f + v , 2 7 +- ^-5 — + VW/ - +-+a —ia —A / = 

sin 77 cos 77 dr\ \ arj J \ sin 77 cos 77 e sin 77 cos 77 4 / 

(3.55) 

Before solving this equation, let us consider the fermions, as we will see their components satisfy a 
very similar equation. 

3.4.2 Fermions 

The operator in this case is 

F =ip-ia (3.56) 
It is convenient to use a left-invariant vielbein. One computes this in toroidal coordinates as: 



d d d 

sin(0 + 4>)^— + cos(6> + 0)(cot?7— - tan 77— ) i = 1 
077 00 0(f) 

d d d 

< - cos(6> + 4>) — + sin(0 + 0) (cot 77— - tan 77— ) i = 2 

OT] 00 Ocp 



Then the Dirac operator can be written as: 



if- 



which, in toroidal coordinates, becomes: 



d_ d_ 
d6 + ~b~4> 

id 3 - § id\ + d 2 
id\ - d 2 -id 3 - I 



i = 3 



(3.57) 



(3.58) 



ip ■ 



.d_ ._a__3 

l de +l dj~ 2 

d .( d 



dr] 



-i(B+<t>) 



+ M cotr ?T^ — tan 77 
06 



d ( d d\\\ 
_ + ^ cotJ7 __ tan?7 _l 

_,d_ 3 
% d0 l d0~2 ) 



(3.59) 

We should also couple to the gauge field, which amounts to the replacement — > + iqG, — > 
-§r — iq (G — 1), as well as the constant a a . Then if we look for eigenspinors of the form 



yj = e 



imO-\-iri(p 



V>i (77) 



(3.60) 



we get the following coupled first-order equations for ipi and ip2'- 
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(4- — (m + qG) cot r/ + (n — <?(G — 1)) tan 77)7/11 = — (m + n + q — A — £cr + 77)7/2 (3.61) 
drj 2 

d 3 

(— (m + qG + 1) cot 77 + (n - g(G - 1) + 1) tan 77)^2 = (m + n + q + \ + ia + -)ipi 

dr] 2 

Solving for ip2 using the first equation and plugging into the second one, we get a second order equation 
in terms of ipi alone: 

— -j- — (m + qG+1) cot r) + (n—q(G— 1) + 1) tan 77 J ( — (m + qG) cot 77 = (n— q(G— 1)) tan 77 ) -01 = 



cf/7 / yd?/ 



(3.62) 



1 3 

= -(m + n + q - A - ia Q + -){m + n + q + A + ia + T^)^ 1 - (3.63) 

This can be rearranged to 

1 d ( #A / (m + g G) 2 (n - g(G - 1)) 2 gg^T/e) ,^2 , ^ , n 

Sin 77 COS ?7— — + — 1 = 1 ; (A+7(7 + -) +1 \tpi = 



0111 1/ wo 1 1 ,11 .9 I 9 I \ ' w Q ' t~\ 

sin 77 cos 77 drj \ drj J \ sin 77 cos z 77 e sm 77 cos 77 z 

(3.64) 

which is precisely the same equation as satisfied by the bosonic modes, with the following relation 
between the eigenvalues: 

A B + (ur + \f = (X F + i°o + \? (3.65) 

This is a quadratic equation for Xf with two solutions Xf±- Their product can be read off as the 
constant term, and one can see that this is simply Xb- 

In principle, we can take i/ii to be the scalar eigenfunction for a given choice of Xb, take one of 
the solutions Xf± to the equation above, and plug these into the first equation in (3.61) to solve for 
ip2- For each Xb and pair of solutions Xf± for which this procedure goes through, we can see that the 
contribution to the partition function, Xf+Xf-/Xb = 1, is trivial. 

However, there are two exceptions we must be more careful with. First, we must also allow 
solutions with ipi = 0, but 7/2 non- vanishing. Then we see that (3.61) reduces to: 

(m + n + q — A — ia a + -)ip2 = 

(-— - (to + qG + 1) cot 77 + (n - q(G - 1) + 1) tan 77)^2 = 
GE77 

Thus the eigenvalue in these cases is X = m + n + q — i cr + ^, and t/> 2 satisfies: 

d , / t \ / .,\ ,\, f —{m + 1) cot 77 + ... near 77 = 

— \0gtp2 = -(m+qG+l) cotr]+(n-q(G-l)+l) t&nr] = < . ^ \ , 1W . , ' „ 

dr] bT y ' ' v v ; y ' [ -(777 + 17+ 1) cot 77+ (77 + 1) tan 77 m the bulk 

(3.66) 

_ J sin~ (m+1) 77 + ... near 77 = . . 

^ V ' 2 ~ 1 sin-( m+9+1 ) 77 cos-^+D 77 in the bulk ^ <] 
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Regularity at the endpoints implies that m and n should be negative integers. But then for —(q + l) < 
m < 0, the bulk solution is singular as it approaches the loop. We will return to this point in a 
moment. These solutions correspond to extra fermionic modes that we have not accounted for before, 
so their eigenvalues should be included in the numerator of the partition function. 

The other exception occurs when the differential operator acting on ipi in the first equation in 
(3.61) annihilates our choice of ipi. Then we must pick ip2 = 0, and there will not be two choices of 
Xf, but only one, and so the cancellation with the corresponding bosonic mode will not be complete. 
We see that in this case (3.61) gives: 

(-^- - (m + qG) cot r] + (n - q(G - 1)) tan 77)^1 = 
dr\ 

3 

(m + n + q + \ + %a + -)^i = 
Now we find A = — (m + n + q + ia + §), and ipi satisfies: 



d i / t , n\ +. i <n + / m cot 77 + ... near 77 = 

— log ib\ = (m + qG) cot 77 — (n — q(G — 111 tan 77 = < . , , „ 3.68 

drj BV v H ' ' K yv " 1 \ (m + q) cot 77- n tan 77 in the bulk v 7 

J sin" 1 77 + ... near 77 = 

\sin( m +^77Cos"77 in the bulk 1 J 

Now regularity at the endpoints forces m and n to be nonnegative integers, and the bulk solutions 
are singular for < m < — q. These are modes for which there is only one solution, say \f+, to 
the equation (3.65), and so the cancellation with the bosons is not complete. The net contribution is 
Af+/Ab = 1/Ajr_, and one can read this off from (3.63) as Xf- = m + n + q — ia Q + ^. 

Putting this together, we see all modes cancel out of the partition function except for the special 
cases noted above, and these give: 

Z = n -<°; m + n + g ~' g °~fi (3.70) 

\Am,n>o( m + n + V - l(T o + 2) 
I) X r, 3 \ t + l 



TT '"° x 2 (3.71) 



(3.72) 



Thus we have obtained the result for the 1-loop determinant without ever having to solve the 
second order differential equation. However, note that this function may have poles, e.g. when q = i. 
To get a better understanding of the origin of these singularities, we will now solve the differential 
equation explicitly. 

3.4.3 Solving the eigenvalue equation 

Consider the following second order equation, which has come up for both the bosons and fermions: 

1 df. df\ ( {m + gGf {n-q{G-l)f qg(r,/e) f ,,^\ f n 

sin 77 cos 77— ] + ( — 1 ^ 1" 7; +1— {£+1) ) / = 



sin 77 cos 77 dr] \ dr\ ) \ sin 2 77 cos 2 77 e sin 77 cos 77 



(3.73) 
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where, for later convenience, we have written the eigenvalue in terms of a parameter t. This is related 
to the eigenvalues by: 



X B = (e+l) 2 -(ia + ^) 2 , \ P ±=±(t + l)-ia -± (3.74) 

Near Loop Region Let us start by focusing on the region < r\ < e, as this will determine what 
boundary conditions to impose on the bulk solution. We start by defining £ = rj/e, and expanding the 
equation above to leading order in e: 

_i^| ) + ( (^ + ?f >^ 

One solution can be immediately obtained, independent of the functional form of g and G, by noting 
that this equation follows from the first order equation: 

| = ^/ (3.76) 

as can be easily checked using G' = g. Then, using G — > as \ — > and G — > 1 as \ ~ ► 00 ; we see 
that: 



3 ~ * ^ £1 



For m > 0, this is the regular solution we want. We can match with the solution in the bulk by looking 
at the behavior in the the large £ region, and we see that we should take the bulk solution which goes 
as f(rj) ~ 7 y( m +9). Note that if q < 0, for < m < —q, the bulk solution would appear to be singular 
right up until we reach the near loop region, at which point the presence of the defect modifies the 
solution to go as r\~ m and be regular. For m > —q, we take the regular solution, as in the absence 
of a defect. Note that this behavior does not depend on the precise functional form of g, only that it 
correctly reproduces a delta function in the e — > limit. 

It remains to understand what happens when m < 0. Here, the first order equation does not have 
any non-trivial solutions regular at £ = 0, so we must return to the second order equation. Then we 
can find a regular solution, but it appears it depends non-trivially on g and G. Specifically, we find: 

f c n £ -> o 

"® ~ \ A(g, m, + B(g, m, g)£"l m +«l £ ^ oo (378) 

for some constants which depend non-trivially on g, m, and q, and in particular are generically both 
non-zero, unlike in the previous case. 

Nevertheless, if we reinstate the e dependence, we see that to match with the bulk, we should take 
the solution there to go as: 

f(ri) - A(g,m,q)^ m+ ^e^ m+ ^ + B(g 7 m 7 q)i-i^ m+ ^^ m+ ^ (3.79) 

and for e — > 0, the first term will dominate, and so we should take the regular solution. 

Thus we should always take the regular solution, / ~ r^l Tn +9l j except in the case where the 
coefficient of the regular part is precisely zero, which happens only when < m < — q, and in 
these cases we take the singular solution / ~ jj _ l m +«l. 
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Bulk Since G = 1 and ,9 = everywhere except an infinitessimal region near r\ = 0, in the bulk the 
equation reduces to: 

- srnr, cos r,f) + ( + + 1 - (* + I) 2 V = (3.80) 

sin 77 cos r\ dr]\ dr\ ) \ sin r\ cos 2 77 / 

It is convenient to look for an / of the form: 

f(r)) = sin™ n cos" V h{sm 2 77) (3.81) 
for m,h we will choose in a moment. Plugging this in, we find the following equation for h(x): 

x{l-x)h" + (m+l-(rh + fi+2)x)h' -\( {m + g) " ~ — + {U - ^ - + (m + n+l) 2 - (i+l) 2 )h = 

(3.82) 

The greatest simplification is achieved by setting rh 2 = (m + q) 2 and n 2 , with the sign of m and n to 
be fixed later. Then we are left with the hypergeometric equation: 

a;(l - x)h" + (c - (a + b + l)x)ti -abh = (3.83) 

where: 

c = 7?i+l, a + b = m + h + l, ab = -((m + n + l) 2 - (£ + l) 2 ) (3.84) 

=>a,6 = i(m + ?i + l±(£+l)) (3.85) 
The solutions can be written in terms of hypergeometric functions 2-^1(0, b; c; x): 

2 F 1 (a J 6;c;x) = f;^%x" (3.86) 

where (a)„ = a(a + l)...(a + n — 1). Provided this is well-defined 2 , it converges for all \x\ < 1, and for 
\x\ = 1 if Re(c-a-fc) > 0. 

In general, there are two linearly independent solutions, but if we impose regularity at the bound- 
ary at x = 1 this restricts us to the solution: 

h(x) = 2 Fi (a, b; a + b + 1 - c; 1 - x) (3.87) 



r(o + 6 + i-c)r(i-c) , , , r(a + b + 1 -c)r(c- 1) , „ „. 

= T^Tl ( 2Fl a ' 6; c; x +- 2^ a + l-c, 6+1-c; 2-c; x 

1 (a + 1 — c)l (0 + 1 — c) 1 (a)l (0) 

(3.88) 

= Ci (a, 6, c) 2-F1 (a, 6; c; x) + C 2 (a, 6, c) 2 -Fi(a + 1 - c, b + 1 - c; 2 - c; x) (3.89) 

we should also set h = \n\ to ensure f{rf) is regular. 

Now we need to fix the behavior at the other endpoint. We have: 



2 Specifically, if c is a negative integer, we must have a or b to be a negative integer greater than or equal to c for this 
series to make sense. 
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/(»7)~Ci(a > 6,c) f7 * + C7a(a,6 > c)»/-* (3.90) 

Thus, depending on the behavior we want here, we should fix the relative values of the Ci which will 
impose a condition on £ and restrict us to a discrete spectrum. Actually, since we are free to choose 
the sign of rh, we can simply pick rh so that the behavior we want is ry Tn , and then the condition is 
always C2 = 0. This is ensured if the argument of one of the Gamma functions in the denominator in 
C2 is set to a negative integer, which determines: 

a,b= -(to + \n\ + l±(£ + 1)) = -k, k e Z >0 (3.91) 

=> X B = {£ + if + K - ^f = (rh+ \n\ +2k + if + {a a - \f (3.92) 

Xp = ±(l + 1) - ia a - i = ±(m + |n| + 2fc + 1) - ia - - (3.93) 
Finally, we recall that the correct choice of rh was found above to be: 



(3.94) 



-\m — q\ < m < —q 
\m — q\ else 

Then the set of eigenvalues is given by taking the expressions above for all integers m, n and non- 
negative integers k. 

Actually, as above, we need to be more careful with the fermions. Specifically, there will be 
eigenvalues in addition to these, as well as some that we have to throw out, corresponding to the cases 
where, respectively, the top and bottom components of the fermion vanish. These are precisely the 
exceptions noted in the previous section, from which the entire contribution to the partition function 
comes. 

One interesting property of these eigenvalues is that, since rh may be negative, the eigenvalues 
may become zero, and the bosonic eigenvalues may even be negative. The first place these zero modes 
can occur is for q = — s, in which case for m = k = a = 0, Xb and are both zero. Actually 
this fermionic eigenvalue is one of the spurious ones that we should throw out, and so in fact there 
is a single bosonic zero mode which results in a pole in the 1-loop partition function. For larger q, 
one finds negative bosonic eigenvalues, which become difficult to make sense of in the path integral. 
We will take the viewpoint that the vortex loop operator is only properly defined for q > — | (and, 
when the matter content is in a self-conjugate representation of the relevant flavor symmetry, this also 
forces q < 5), although the naive result for the determinant gives a natural analytic continuation of 
this quantity which may have some physical relevance. 



4 Duality with Vortex Loop Operators 

Let us now turn to some applications of the vortex loop operator. Since this operator is defined by an 
explicit procedure applied to a global U(l) symmetry, then if we know what this symmetry maps to 
across a duality, we obtain an identification of loop operators on each side of the duality. This provides 
a new operator mapping across duality, although it does not provide any essentially new information 
beyond the mapping of global symmetries. 

One inter?sting case is when a U(l)j symmetry on one side of the duality is identified with a flavor 
symmetry on the dual. This occurs, for example, in 3D mirror symmetry. Let us take the simplest 
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case, the duality between M = 4 SQED with one flavor and a free hypermultiplet. Note that there 
are Wilson loops on the SQED side, but it is less clear what the corresponding loop operators on the 
dual side are, since there is no gauge group. 

However, the results above provide the answer: the Wilson loop on the SQED side are the same as 
defects in the global U(l)j symmetry. Mirror symmetry dictates that this symmetry is identified with 
a £7(l)y flavor symmetry under which the two chirals in the hypermultiplet have charge ±1. Then 
the Wilson loop in SQED simply maps to a defect operator in this flavor symmetry, supported on the 
same loop. 

At the level of the matrix model, this follows from the identity: 

/ 2-n-ir/A p 2irq\ i 
d\- — fr— = - — — (4.1) 
2 cosh(7r A) 2 cosh 7r(?y — iq) 

which simply follows from extending the usual self- Fourier-transform property of 1/ cosh to the entire 

complex plane. 

It is worth noting that the integral on the LHS only converges for |g| < i. As remarked above, 
this is precisely the range in which the defect operator is also well-defined. The divergence on the LHS 
as |g| — > \ is reflected on the RHS as a bosonic zero mode developing in the defect background. The 
RHS gives a natural analytic continuation of this quantity to \q\ > i, but it is not clear what, if any, 
physical relevance this has. 

We can also consider mirror symmetry applied to SQED with Nf flavors. Here the dual is a quiver 
theory with gauge group U(l) Nf /U(l)di ag — £7(l) iV -'' _1 , with Nf bifundamental flavors (q a q a ) charged 
as (1,-1,0, ...), (0, 1, —1, ...), (—1, 0, 1), Nf neutral chirals S a , and a superpotential ^2 a q a S a q a - 
Here the £7(1) j symmetry of SQED maps to a flavor symmetry under which the all the q a have charge 
■j^-. In particular, the Wilson loop in SQED maps to a defect in this symmetry. Note that, because of 
the extra factor of jfj in this mapping, the restriction on the defect charge now allows us to consider 

Wilson loops of charge up to ±-^. This coincides with the fact that, in the matrix model, the Wilson 
loop expectation value now converges for this wider range of charges because of the increased damping 
in the integrand from the factor (2 cosh(7rA)) _JV/ . 

Finally, we note that we can also apply this type of argument to theories with U(N) gauge 
symmetries for N > 1, provided we restrict only to Wilson loops in the overall £7(1) of the gauge 
group. For example, in the U(N) versions of Aharony ([18]) and Giveon-Kutasov dualities ([19]), it is 
known that the £7(1) j current maps to itself, up to a flip of sign, and so a Wilson loop in the overall 
£7(1) of the gauge group must map to the same Wilson loop, with the opposite charge. One can also 
consider such abelian Wilson loops in non-abelian mirror symmetry, where they will map to flavor 
defects as above. 

5 Discussion 

We have defined a set of abelian vortex loop operators which exist in any abelian gauge theory in 2 + 1 
dimensions. The definition can be extended to conformal field theories with abelian global symmetries 
by weakly gauging the symmetry currents. Witten's SL(2, Z) action for this class of theories naturally 
extends to loop operators. In fact, abelian vortex loops are the S-duals of the ordinary Wilson loops. 
Alternatively, a Wilson loop in a £7(1) factor of the gauge group can be viewed as a global vortex loop 
for the associated topological £7(1) j symmetry. 

One of our results was the definition of the supersymmetric version of the vortex loop. This was 
accomplished by identifying the fields in the abelian vector multiplet which needed to be turned on 
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to create the right type of singularity. It turns out that, besides the singular gauge connection, we 
needed to also give an imaginary background value to the auxiliary scalar D. This is an interesting 
example of the fact that background fields need not satisfy the reality conditions usually imposed on 
the dynamical fields of the theory We proceeded to evaluate the expectation value of a supersymmetric 
defect loop, defined on a great circle on S 3 , using localization. The result could be anticipated by 
considering the SX(2,Z) action and, indeed, had already been derived from the 4d perspective. We 
have given, by using and comparing two different regularization methods, an additional microscopic 
derivation. 

The supersymmetric vortex loop plays a central role in mirror symmetry of 3d gauge theories. This 
class of dualities has the property that it exchanges flavor symmetries with the topological symmetries 
associated to the abelian factors of the dual gauge group. As a consequence, the duality exchanges (the 
supersymmetric versions of) U(l) Wilson loops with abelian vortex loops. Identifying such entries in 
the duality dictionary is an important step towards, possibly, proving the duality for the full quantum 
theory. We have demonstrated that the expectation values for the dual loop operators match, in simple 
examples, by using localization and the matrix model. 

The analysis presented here has a natural extension to non-abelian defects. The definition of 
such an operator can require additional steps to ensure gauge invariance. For a defect in a global 
non-abelian symmetry group the definition is similar to the abelian case and the results can be read 
off from Section 3 by conjugating the defect data into the Cartan of the flavor group. When the defect 
appears in a dynamical gauge field, the localization procedure for the vector multiplet is modified. The 
result in the case of pure Chern-Simons theory is known to require a quantization of the data entering 
the definition of the defect, the overall effect being, again, the insertion of a Wilson loop operator in 
some representation [6]. A naive analysis would imply that this result is not affected by the presence 
of additional charged matter. However, the mapping of such operators under mirror symmetry for 
non-abelian theories and under Seiberg-like dualities requires further investigation. We will pursue 
this investigation elsewhere [? ]. 
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